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Abstract — A geometric form of Euler-Lagrange equations 
is developed for a chain pendulum, a serial connection of n 
rigid links connected by spherical joints, that is attached to a 
rigid cart. The cart can translate in a horizontal plane acted 
on by a horizontal control force while the chain pendulum 
can undergo complex motion in 3D due to gravity. The 
configuration of the system is in (S 2 ) n xR 2 . We examine the rich 
structure of the uncontrolled system dynamics: the equilibria 
of the system correspond to any one of 2 n different chain 
pendulum configurations and any cart location. A linearization 
about each equilibrium, and the corresponding controllability 
criterion is provided. We also show that any equilibrium can 
be asymptotically stabilized by using a proportional-derivative 
type controller, and we provide a few numerical examples. 

I. Introduction 

Pendulum models have been a rich source of examples in 
nonlinear dynamics and control [1], [2]. For example, the 
dynamics of a double spherical pendulum and a Lagrange 
top have been studied in [3], [4], [5]. Generalized models, 
such as a 3D pendulum [6] or a 3D pendulum attached 
to an elastic string [7], have been considered. A variety 
of control techniques have been applied, such as passivity- 
based approaches [8], [9], swing-up strategies [10], [11], 
Lyapunov-based method [12], controlled-Lagrangian [13], 
and hybrid control systems [14]. In particular, stabilization 
of a triple inverted pendulum has been studied in [15], [16]. 

In this paper, we consider the dynamics and control of a 
chain pendulum on a cart, that is a serial connection of n 
rigid links, connected by a spherical joint, attached to a cart 
that moves on a horizontal plane. The configuration manifold 
of this system is (S 2 )™ x R 2 , where the manifold of unit 
vectors in M 3 is denoted by the two-sphere S 2 . 

Many interesting mechanical systems, such as robotic 
manipulators or variations of pendulum models, evolve on 
the two-sphere S 2 or on products of two-spheres (S 2 )". In 
most of the literature that treats dynamic systems on (S 2 ) n , 
including many of the above references, either 2n spherical 
polar angles or n explicit equality constraints that enforce 
unit lengths are used to describe the configuration of the 
system. These descriptions necessarily involve complicated 
trigonometric expressions and introduce additional complex- 
ity in analysis and computations, as well as singularities. 
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Fig. 1. Chain pendulum on a cart (n = 3) 

We demonstrate that globally valid Euler-Lagrange equa- 
tions on (S 2 )™ x M. 2 can be developed for the chain pendulum 
on a cart system, and they can be analyzed in a compact form 
without local parameterization or constraints. This also leads 
a coordinate-free form of linearized equations, controllability 
criteria, and control systems. The main contribution of this 
paper is providing an intrinsic and unified framework to 
study dynamics and control of chain pendulum on a cart 
systems, that is uniformly applicable for an arbitrary number 
of links, and globally valid for any configuration of the links. 

II. Lagrangian Dynamics of the Chain Pendulum 
on a Cart System 

A. Background 

The cart of mass m can translate on a horizontal plane, 
and its position in an inertial frame is denoted by i € R 2 . 
It is acted on by a horizontal control force net 2 . A serial 
connection of n rigid links, connected by spherical joints, is 
attached to the cart, where the mass of the i-th link is denoted 
by nii and the link length is denoted by For simplicity, 
we assume that the mass of each link is concentrated at the 
outboard end of the link. The direction vector of each link in 
an inertial frame is given by qi e S 2 = {q G M 3 | ||g|| = 1} 
for i = 1, . . . , n. The configuration manifold of this chain 
pendulum on a cart system is (S 2 )™ x K 2 . 

The inertial frame is defined by the unit vectors e\ = 
[1; 0; 0] G K 3 , e 2 = [0; 1; 0] e M 3 , e 3 = [0; 0; 1] G M 3 ; 
we assume that e 3 is in the direction of gravity. Define C — 
[ei, e 2 ] g M 3x2 . 

The kinematic equation for the direction vector of the z-th 
link qi is given by 

qi = Ui x qi = Uiqi, (1) 

where coj G M 3 is the angular velocity of the i-th link 
satisfying uii ■ qi = 0. 



The hat map A : R 3 — > so (3) transforms a vector in R 3 to 
a 3 x 3 skew-symmetric matrix, and is uniquely defined by 
the property that xy — x x y for any x, y £ R 3 . The inverse 
of the hat map is denoted by the vee map V : so (3) — > R 3 . 

Throughout this paper, the dot product of two vectors is 
denoted by x ■ y — x T y for any x, y £ R n . The nxn identity 
matrix is denoted by /„. The ra x to by matrix composed 
of zero elements is denoted by O n xm, and it is written as 
n if n = ra. A column-wise stack of matrices is written 
as [A; B] = [A T 1 B T ] T £ R( tt+i, > xn for A G R ax " and 
B £ R &xri . Some properties of the hat map are given by 



xy = x x y = —y x x = —yx, 
x ■ yz = y ■ zx = z ■ xy, 
xyz = (x ■ z)y - (x ■ y)z, 
C T C 



-C 1 e z 3 C = I 2 



(2) 
(3) 
(4) 
(5) 



for any x,y,z £ R 3 . 
B. Lagrangian 

The location of the cart is given by Cx £ R 3 in the inertial 
frame. Let x$ £ R 3 be the position of the outboard end of 
the i-th link in the inertial frame. It can be written as 



Xi — Cx + ^ ' laQa- 



(6) 



The total kinetic energy is composed of the kinetic energy 
of the cart and the kinetic energy of each mass: 



T = ~rn||x|| 



1 " 



\Cx- 



^2l a qa\\ 2 - 

a=l 



(7) 



For simplicity, we first consider the part of the kinetic 
energy, namely T x that is dependent on the motion of the 
cart. The part of |7} dependent on x is given by 



1 ™ 



n n 



T t = -(to+> m,i) \\x\\ + Cx ■ y y m a liqi 

i— 1 i— 1 a—i 

This can be written as 



1 

T x = -M 00 ||a; 



| 2 +X-J2 M 0r& 

i=l 



where the inertia matrices M 0Q £ R, M oi £ 
M i0 £ R 3x2 are given by 

n n 
i—1 a—i 

(9) 

for i = 1, . . . , n. The part of the kinetic energy, namely T q , 
that is independent of x is given by 



i 



i 

i—1 a—1 *)J = 1 

where the inertia constants Mij € M are given by 



IE- 

a=l 



M,. 




(11) 



for i,j = 1, 



From (|8]l and ( [10) , the total kinetic energy is given by 

Mnnll.xll^ 4- T, ■ 



n 1 n 

T = -M o\\x\\'' M oA + 2 E ' (12) 

where the inertia matrices are given by (|9j, ( [TT} . The grav- 
itational potential energy is composed of the gravitational 
potential energy of each mass. From |6]), it can be written as 

n n n 

V = - ^ m i9Xi ■ £3 = - E E m a9he 3 ■ q t . (13) 

i—1 i—1 a—i 

The Lagrangian of a chain pendulum on a cart is L = T— V . 

C. Eider -Lagrange equations 

In [17], it is shown that a coordinate-free form of Euler- 
Lagrange equations on the two-spheres can be derived from 
Hamilton's variational principle. The key idea is to express 
the variation of a curve on the two-sphere in terms of the 
exponential map on SO(3). Let qi(t) be a curve on S 2 . Its 
variation can be written as 

qt{t)=exp{eZ i {t))q i (t), 

where £j(t) is a curve in R 3 satisfying q^t) ■ ^(t) = 
for all t. As the exponential map represents the rotation of 
qi{t) about the axis by the angle |e|||£i(i)ll at each t, 
it is guaranteed that the varied curve g| (t) lies in S 2 . The 
corresponding infinitesimal variation is given by 

exp(e&(i))«?i(t) = x qi (t), (14) 



Sq z (t) 



dc 



which lies in the tangent space T g .( t )S 2 as it is perpendicular 
to qi(t) at each t. Using these, we obtain the Euler-Lagrange 
equations of a chain pendulum on a cart as follows. 

Proposition 1: Consider a chain pendulum on a cart, 



whose Lagrangian is given by (12i and (13i. The Euler- 
Lagrange equation on (S 2 )" x R <ire as follows: 
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-\\q n \\ 2 M nn q n - m n ^„9 2 e 3 
where the inertia matrices are given by d9) and ( 1 1 



Or equivalently, it can be written in terms of the angular 
velocities as 
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+ YZ=2 m a,ghq2e 3 



qi =w,x q t . 
Proof: See Appendix [A] 



(16) 



(17) 



D. Remarks about the Chain Pendulum on a Cart System 

The system of equations above for the chain pendulum 
on a cart system remains valid for arbitrary deformations 
of the chain with respect to the cart; part or all of the chain 
pendulum may be above the cart while part or all of the chain 
pendulum may be below the cart. We ignore all possible 
collisions of the chain pendulum with itself and of the chain 
pendulum with the cart. 

The chain pendulum with the cart inertially fixed is a 
special case of the model above. If there is a single link, n = 
1, the spherical pendulum is obtained; if there are two serial 
links, n = 2, the double spherical pendulum is obtained. The 
resulting dynamics of the chain pendulum for any number 
of links can be very complex. The Euler-Lagrange equations 
for this case are easily obtained as special cases of the prior 
equations. The cart is included in the system primarily as a 
means of achieving control actuation. 

The planar chain pendulum on a cart is also a special case 
of the above model. Let v G R 3 be a vector normal to e 3 , 
i.e., v ■ e 3 = 0. Suppose that the initial condition is chosen 
such that Wi(0) x« = and <7i(0) • v = for i = 1, . . . , n 
and the control input is given such that v ■ Cu — 0. From 



(16 1, we can show that the motion of links with respect to 
the cart remains confined to the plane normal to v. 

In short, equations (15 1, or (16i and ( [T7] i, provide the 
Euler-Lagrange equations for a chain pendulum or a chain 
pendulum on a cart with an arbitrary number of links. As 
they are developed in a compact, coordinate-free fashion, 
singularities that arise when using local coordinates are 
completely avoided. 

III. Dynamic Properties of the Uncontrolled 
Chain Pendulum on a Cart System 

A. Equilibrium Configurations 

We can determine the equilibria of the uncontrolled chain 
pendulum on a cart system using ([16]) when u = 0. Assuming 
the configuration is constant, the equilibrium configurations 



in (S 



2\n 



are given by the cart in a fixed location and 



the chain pendulum aligned vertically, i.e., 

qi x e 3 = 0, for i = 1, . . . , n. 

Consequently, assuming that we are interested in the 
equilibria where x — 0, there are 2™ possible equilibrium 
configurations that correspond to the case where all n links 
are vertical: % = ±e 3 for i = 1, . . . , n. The equilibrium for 
which all links are aligned with the gravity direction, % = e 3 
for i — 1. . . . , n, is referred to as the hanging equilibrium; 
the equilibrium for which all links are aligned opposite to 



the gravity direction q t = — e 3 for i = 1, . . . , n, is referred to 
as the inverted equilibrium. Two other interesting equilibria 
correspond to the case that all adjacent links point in opposite 
directions, that is • qi + i = — 1. These two equilibria, one 
with the first link pointing in the direction of e 3 and the 
other with the first link pointing in the direction of — e 3 , are 
referred to as folded equilibria since the chain pendulum is 
completely folded in each case. 

We introduce binary variables Sj to describe a specific 
equilibrium. It is defined as 



(18) 



1 if qi = e 3 (along gravity), 

Si = < 

1—1 if qi = — e 3 (opposite to gravity) 

for i = 1, . . . , n. Then, an ordered n-tuple s = 
(si, S2, ■ ■ ■ , s n ), specifies a particular equilibrium. 

B. Linearized Equations of Motion 

The local dynamics near each equilibrium can be studied 
using linearization methods. From ( fl4) , the variation of q^ 
from any equilibrium configuration can be written as 

% = d x e 3 or % = & x -e 3 , 

where £j g R 3 with £j • e 3 = 0. The variation of Ui is 
given by 5ui E M 3 with • e 3 = 0. Therefore, the third 
components of £j and Sui for any equilibrium configuration 
are zero, and they are omitted in the following linearized 
equation, i.e., the state vector of the linearized equation is 
composed of C T ^ € K 2 . As a result, the dimension of the 
corresponding unconstrained linearized equation is 2n + 2. 

Proposition 2: Consider an equilibrium of a chain pendu- 
lum on a cart, specified by s = (si, . . . , s n ) and x = 0. The 
linearized equation of the Euler-Lagrange equation ( [TB"] ) is 



Mx + Gx = Bu, 



(19) 



or equivalently 
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where the corresponding sub-matrices are defined as 

x 9 = [C7 T £i; ... ;C T £„], 
M xx = Moo/a, 

M xq = [~s 1 M 01 e 3 C -s 2 Mo2i 3 C ■ ■ ■ -s n M 0n e 3 C] 
M„* = 
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a=l 

where = SiSj for i,j = 1, . 
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Proof: See Appendix |B| ■ 
The local eigenstructure near each equilibrium can be 
determined from the linearized dynamics. The eigenvalues 



of {19]) are the roots of det[A 2 M + G] = 0. Note that there 
are zero eigenvalues since the first two columns and rows of 
G are zero. These correspond to the cart dynamics. 

For the hanging equilibrium, given by Sj = 1 for all 
i = 1, , . , ,n, the matrix G qq becomes positive-definite. 
This implies that A 2 < always. Then, the stability of the 
nonlinear dynamics ( [To*} is inconclusive from the linearized 
equation (19i as Re [A] = 0. But, it can be shown that the 



hanging equilibrium is stable in the sense of Lyapunov mod- 
ulo the cart position, using the total energy as a Lyapunov 
function: the hanging equilibrium is the local minimum of 
the total energy, which is conserved along the solution of the 
uncontrolled system. For each of the other 2™ — 1 equilibrium 
configurations, there exist a pair of positive and negative 
eigenvalues, which implies that it is an unstable saddle. Het- 
eroclinic connections between the various equilibria provide 
a non-local characterization of the dynamic flow. 

IV. Control Analysis 

Various control problems can be posed for the chain pen- 
dulum on a cart system. For example, feedback control might 
be used to achieve asymptotic stabilization of any of the 
natural equilibrium solutions. From the linearized equation 



(19 1, we find a criterion for controllability as follows. 

Proposition 3: Consider an equilibrium of a chain pendu- 
lum on a cart, specified by s — (sj., . . . , s n ) and x = 0. 
Suppose that the inertia matrix M of the linearized equation 



( 19 1 is positive-definite. Then, the equilibrium is controllable 
if and only if the following subsystem is controllable 



M qq X.g 



G qq = M xq U 



for a control input u e I 2 , or equivalently 

rank[A 2 M„ + G qq , M qx ] = 2n 



(20) 



(21) 



for any A e C. 

Proof: See Appendix [C] ■ 
This proposition states that the controllability of (19i is 
equivalent to the controllability of a reduced system given 



by (20 1. It represents the dynamics of the chain pendulum, 
without the cart, where the input matrix is given by M g2 . 
that corresponds to the inertia coupling between the chain 
dynamics and the cart dynamics. 

If the linearized equation about a specific equilibrium 
configuration is controllable, then we can design a control 
system to asymptotically stabilize that equilibrium. 

Proposition 4: Consider an equilibrium of a chain pen- 
dulum on a cart, specified by s = (s 1; ...,s„) and x = 
0. Assume that the corresponding linearized equation is 
controllable (as characterized by Proposition |3j. The control 
force u is chosen as follows: 

n 

u = -K x x - K x x - ^{X ? .C T (sie 3 x q t ) + K u .C T Ui}, 

i=l 

(22) 



for controller gains K x , K±, K qi , K Ui E 



p2x2 



for i 
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Fig. 2. Uncontrolled response: perturbation from a folded equilibrium 



such that the equilibrium of the controlled system is locally 
asymptotically stable. 

Proof: See Appendix [P] ■ 
This control system can be used uniformly to stabilize any of 
2" equilibrium configurations of a chain pendulum on a cart. 
But, as it is based on the linearized dynamics, the region of 
attraction could be limited. 

V. Numerical Examples 

In the subsequent numerical simulations, we consider a 
chain pendulum model with five identical links, i.e. n = 5. 
This system has twelve degrees of freedom with two force 
control inputs on the cart. The properties of the cart and the 
links are chosen as 



0.5 kg, 



0.1kg, k = 0.1m for % = 1, 



Throughout this section, the following units are used: 
kg, m, sec and rad, unless otherwise specified. 

First, simulation results for the uncontrolled chain pendu- 
lum on a cart system are presented. The initial condition is a 
small perturbation from one of the completely folded equi- 
libria, given by s = (—1,1,-1,1,-1). More specifically, 



x(0) = [0.2; 0.1], x(0) = [0;-0.1], 
t?»(0) = s 4 e 3 = (-l) l e 3 for i = 1,. . .,4, 
q 5 (0) = [sinl°;0;-cosl°], 
Wj(0) = 3x i for i = l,..., 5, 



(23) 



where the fifth link is perturbed by 1° from the equilibrium. 
The corresponding simulation results are shown at Figure 
[2] The given initial condition guarantees that the relative 
motion of the links with respect to the cart always lies in 



l,...,n. Then, there exist values of the controller gains the eie 3 plane, which is depicted in Figure 2(c) Figure 2(d) 
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Fig. 3. Controlled response: asymptotic stabilization of the hanging 
equilibrium s = {1, 1, 1, 1, 1} 



demonstrates the complex transfer of potential energy and 
kinetic energy between the links and the cart as the chain 
pendulum on a cart dynamics evolve. 

Second, simulation results are presented that show the 
response of the chain pendulum on a cart system to a 
feedback control (j22j) that stabilizes the hanging equilibrium 
s = (1,1, 1,1,1). The initial conditions are same as (23 1, and 



the controller gains are chosen from a linear quadratic regula- 
tor with the weighting matrices Q = diag[8/2, hn, 8/2, Iin] 
and R = I 2 - We define the following variables that measure 
the direction errors and the angular velocity errors of links: 



Vii 



n 

Ei 



Figure [3] illustrates that the chain pendulum on a cart 
asymptotically approaches the hanging equilibrium. 



Next, we consider the control system ( 22 1 to stabilize a 
partially folded equilibrium given by s = ( — 1,-1,-1,1,1), 
i.e., the first three links are opposite to gravity, and the 
remaining last two links are aligned with gravity. Initial 
conditions are given as follows: 

gx (Q) = -[sin6 o ;0;cos6°], 
9 2 (0) =<fe(0) = -[0;sin4 o ;0;cos4°], 
q 4 (0) = [sin 5° cos 4°; -sin 5° sin 4°; cos 5°], 
q 5 (0) = [-sin35°;0;cos35°]. 

Other initial conditions for x(0),x(0) and Wj(0) are iden- 



tical to (23 1. Controller gains are chosen from a lin- 
ear quadratic regulator with the weighting matrices Q — 
diag[/ 2 , Shn, Hi 8/2n] an d R — h- Figure [4] illustrates that 





- 20 o L 

40 r 
20 




2 4 6 8 10 

t 

(d) Control force u 



(c) Location of with respect 
to the cart (J2i=i 



Fig. 4. Controlled response: asymptotic stabilization of a partially-folded 
equilibrium s = { — 1, —1, —1, 1, 1} 



the cart and the pendulum asymptotically converge to the 
equilibrium s = ( — 1, —1, —1, 1, 1). 

VI. Conclusions 

Euler-Lagrange equations that evolve on (S 2 )" x M 2 have 
been derived for the chain pendulum on a cart system. These 
equations of motion provide a remarkably compact form 
of the equations of motion, which enables us to analyze 
their dynamic properties and control systems uniformly for 
an arbitrary number of the links, and globally for any 
configuration of the links. 

We emphasize that modeling, analysis, and computations 
can be carried out directly in terms of a geometric coordinate- 
free framework as illustrated by the chain pendulum on a 
cart system studied in this paper; there is no need ever to use 
local angle coordinates. This important fact is not appreciated 
by many researchers in dynamics and control who continue 
to formulate many dynamics and control problems in local 
coordinates, which could otherwise be analyzed with greater 
ease in a coordinate-free framework. 

Appendix 

A. Proof of Proposition [7] 

The variations of the Lagrangian with respect to x and x 
are given by 

T> X L ■ 5x = 0, 

n 

T> ± L ■ Sx = (M QO x + 2J M 0i qi) ■ Sx, 
»=i 

where T> X L represents the derivative of L with respect to x. 
From (14 1, the variation of qt is Sqt = £i X qi for G M 3 



with £j • qi — 0. The variation of the Lagrangian with respect 
to qi is given by 

n n 

T> qi L ■ Sq t = '^2 m a gke 3 ■ (& x q, t ) = - } m a gke 3 qi ■ 

a— i a—i 

where Q has been used. The variation of qi is given by 

S<}i =& X qi +£ x <ii- 

From this and ([3]), the variation of the Lagrangian with 
respect to <ji is given by 

n 

T> 4i L ■ 5q\ = (M i0 ± + ^2 M ijQj) ■ (ii x 1 + & x Qi) 

3=1 

n n 

= qi(M lQ x + ^ Mijqj) ■ ii + q t (M lQ x + ^ M ijQj) 1 £«• 
i=i i=i 
Using these expressions and integrating by parts, the 
variation of the action integral can be written as 

8<8= / -{M ox + 22M 0l qi}-5x 

n n n 

+ y^{-gi(Mj x + Mij^,-) - m a glie 3 qi} ■ & dt. 

z—l J— 1 a—i 

According to the Lagrange-d'Alembert principle, the sum 
of the variation of the action integral, and the integral of the 
virtual work done by the control force on the cart, namely 
J t f u5x dt, is zero. This implies that the expression within 
the first pair of braces in the above equation is equal to 
—u, and the expression within the second pair of braces is 
parallel to for any 1 < i < n, as is perpendicular to q^. 
Therefore, we obtain 



(24) 



(M l0 x + M ijQj) + J2 m *9hq>e 3 = 0. (25) 

3=1 o=t 



Equation ( 25 1 is rewritten to obtain an explicit expression 
for q\. As qi - qi — 0, we have qi ■ <ji + qi - q\ — 0. Using this 
and Q, we have 

-QiQi = "{ii ■ Qi)Qi + {Qi ■ Qi)Qi = (?« ' Qi)Qi + Qi- 

Substituting this into p5j), 



These equations ( |24) i and ( |26] i are rewritten in a matrix form 
to obtain ( fT5] l. 

These can also be rewritten in terms of the angular 
velocities. Since q\ = u>i x qi for the angular velocity 
satisfying qi ■ u)i = 0, we have 

q\ = w, x qi + iOi x (uj, x q t ) = -q i uj l - \\oj l \\ 2 q i . 



MuQi - Qii M ioZ + J2 M ijQj) 

3=1 

n 

= -Mu\\qi\\ 2 qi - ^2m a gkq 2 e 3 . (26) 



Substituting this into ( |24|) and ( |26) i, and using the fact that 
• Qi = 0, we obtain ( |16| l. 

B. Proof of Proposition [2] 

Consider the hanging equilibrium where s = (l,l,...,l), 
and x = 0. The variations from the hanging equilibrium are 

x e = eSx, i £ = eSx, q\ = exp(e^)e3, u\ = eSuJi, 

where Sx,Sx € M 2 , and ^i,Suji g ]R 3 with & • e3 = and 
(5wi • e3 = 0. This yields the following infinitesimal variation 
<5<7i = £j x e 3- From ([TJi, 5gj is given by 

&Qi = C, x e 3 = x e 3 + x x e 3 ) = 5a;, x e 3 . 

Since both sides of the above equation is perpendicular to 
e 3 , this is equivalent to 63 x x 63) — e 3 x (6u!i x 63), 
which yields 

£ - (e 3 • £)e 3 = 5wj - (e 3 ■ 5w,)e 3 . 

Since ^ • e 3 = 0, we have £ ■ e 3 = 0. As e 3 ■ SuJi = from 
the constraint, we obtain the linearized equation for ([TJ: 



£,t = Suji 



(27) 



Substituting these into ( [To} , and ignoring the higher order 
terms, we obtain 



-^00-^2x2 

C T e 3 M w 
C T e 3 M- 



211 



-M 01 e 3 C 
M 2 ih 



C 1 e 3 M n0 M nl I 2 



-M 02 e 3 C 
M 12 I 2 
M 22 I 2 

M n2 I 2 



-M 0n e 3 C 
M ln I 2 
Manh 





Sx 




u 










- Eo=i m aghd 




X 






~ Ea=2 m agH' 2 


(28) 




c T 'L. 




_ -m n gl n £' n _ 





where we have used the fact that e 3 — diag[— 1, — 1, 0], 
C T e 3 C = —I 2 and e 3 CC T — e 3 . This is the linearized 
equation about the hanging equilibrium. 

This analysis can be easily generalized to other equilibria, 
where one or more links are aligned opposite to the gravity. 
Consider the equilibrium where only the z-th link is pointing 
upward, i.e. tfo = e 3 and qj — —e 3 for all j ^ i. By 
following the same procedure, we obtain the same form of 
the linearized equation as ( |28] i, where all of the terms related 
to Mij, Mji and /j for all j 7^ i are multiplied by —1. This 



yields (19 1 



C. Proof of Proposition [i] 

Suppose that M is invertible. It is well-known that the 
linearized system ( fT~9| > is controllable, if and only if 

rank [A 2 M + G, B] = In + 2 (29) 

for any generalized eigenvalue A satisfying det[A 2 M + G] = 
(see [18], [19]). This is a generalization of the Popov- 
Belevitch-Hautus (PBH) eigenvalue test to a second-order 



system. While it is not explicitly stated in the above refer- 
ences [18], [19], it is straightforward to find an equivalent 
condition in terms of eigenvectors, which is similar to the 
PBH eigenvector test. 



where K x 

\K X ,K U . 



...,K t 

i2x2n+2 



<l„l 



e 



p2x2n+2 



and K-o 



m , . . . , n Wii j c ia. . Since ( 19 1 is controllable, 

we can choose the controller gains K x 



, Ky such that the 



We claim that (29) holds if and only if there is no general- 
ized left eigenvector that is orthogonal to B, i.e. for any non- 



equilibrium is asymptotically stable for the linearized equa- 
tion (fl9j. According to Theorem 4.7 in [20], the equilib- 
rium becomes asymptotically stable for the nonlinear Euler- 



zero eigenvector Vj G 



p2n+2 



satisfying v/ (A 2 M + G) = 0, Lagrange equation ( 16 1 



we have vf B ^ 0ix2- The proof is as follows: 
(Sufficiency) Suppose that there is a generalized eigenvector 



Vj that is orthogonal to B. Left-multiplying (29i by Vj yields 

vf [A 2 M + G, B] = [vf (A 2 M + G), vf B], 

which becomes [0i x2 n+2, 0i X 2] when A = Aj. Therefore, 
the matrix given in p9] > has linearly dependent rows, which 
implies that it is rank-deficient. 



(Necessity) If the matrix given in (29 1 is rank-deficient for 
Aj, there exists a vector Vj satisfying 

vf [A 2 M + G, B] = [((A 2 M + G)v t ) T , vf B] = [0 lx2n+4 ], 

which implies that Vj is a generalized eigenvector that is 
orthogonal to B. 



Using this eigenvector test, we show that (21 1 implies d29 



More specifically, we show that if ( 29 1 is false, then pTj) 
is false. Suppose that there exists a generalized eigenvector 
Vj = [v x ; v g ] that is orthogonal to B. Then, we have vf B = 
r£ — 0i X 2 from the definition of B in 



v f 2nX 2 



(19 1. As Vj is the left eigenvector, we also have 

A?M, 



v/[AfM + G] = [0^ 1 ,v ] 



[A 2 v^M 



[0 



1x2 



qx * q 
0lx2n • 



(A 2 M 9 



A 2 M 9:c 
+ G, 



A?M 



qq 



G 



</'/. 



<i'i> 



(30) 



= 2 



When A, = 0, this yields vf G qq = 0i x2 „ =5 
as G qq is invertible. This is not possible since it contradicts 



the fact that v t = [0 2x i;v 9 ] ^ 2n+2 > 
Then, ([30| implies 



Therefore, Aj ^ 0. 



V^ (XfM qq + G qq ) = lx2 n, V q M qx =0 lx2 , (31) 

which states that there exists a non-zero generalized left 
eigenvector of (20i that is orthogonal to M. qx . Therefore (21 1 
is false. 

As a last step, we show p9} implies pi) . If pT) is false, 
there exists a non-zero eigenvector v 9 satisfying pT) . Let 
Vj = [02xii v g ]. Then, it is orthogonal to B as vf B = 0i x2 - 
And Vj is a generalized left eigenvector of (M, G) as it 
satisfies (30i. In short, pi) is equivalent to p9| . 



D. Proof of Proposition |4] 



Using ( 14 1, ( 27 1, the linearized control input is given by 

n 

u = -K x Sx - K x Sx - y^{K qt C T £i + K Ui C T 5u>i}. 



=1 



From the 



definition of the state vector x = 
■ ; C T Cn] m ( 19 1, u can be written as 



u = K x x — KxX, 
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